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ABSTRACT 


Gallium  Nitride  (GaN)  High  Electron  Mobility  Transistors  (HEMT’s)  are 
microwave  power  devices  that  have  the  performance  characteristics  to  improve  the 
capabilities  of  current  and  future  Navy  radar  and  communication  systems.  The  Office  of 
Naval  Research  (ONR)  is  funding  research  for  the  development  of  GaN -based 
microwave  power  amplifiers  for  use  in  future  radar  and  communication  systems.  This 
thesis  studies  the  effects  of  AlGaN/GaN  HEMTs'  polarization,  piezoelectric  (PZ)  and 
spontaneous,  properties  utilizing  the  commercially  available  Silvaco  Atlas™  software  for 
modeling  and  simulation.  The  polarization  properties  are  suspected  to  enhance  the  two- 
dimensional  electron  gas  (2 DEG)  at  the  AlGaN/GaN  interface. 
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EXECUTIVE  SUMMARY 


A  computer  model  of  Aluminum  Gallium  Nitride  (AlGaN)/Gallium  Nitride 
(GaN)  High  Electron  Mobility  Transistors  (HEMTs)  has  been  developed  for  the  purpose 
of  designing  and  modifying  high  power  and  high  frequency  power  amplifiers.  GaN- 
based  HEMT  power  amplifiers  are  being  investigated  for  application  in  future  military 
radar  and  communication  systems.  This  model  includes  the  impact  of  the  AlGaN/GaN 
polarization,  piezoelectric  and  spontaneous,  effects  based  upon  the  material  parameters  of 
AlGaN,  GaN,  and  HEMTs. 

The  polarization  effects  result  from  the  GaN's  noncentrosymmetric  crystal 
structure  combined  with  strain  induced  at  the  AlGaN/GaN  heterojunction  due  to  their 
lattice  mismatch.  This  effect  results  in  a  large  electron  concentration  at  the 
hctcroj unction  thus  providing  exceptional  conductivity  of  electron  current.  As  a  material, 
GaN  has  a  large  energy  bandgap,  high  breakdown  voltage,  and  high  peak  electron 
velocity.  The  HEMT’s  characteristic  heterojunction  produces  a  quantum  well  which 
further  enhances  the  speed  of  electron  transport  from  source  to  drain  contacts  and  reduces 
electron  scattering. 

The  Silvaco  software,  a  physics-based  modeling  program,  was  utilized  to  model 
and  simulate  an  AlojGao.vN/GaN  HEMT.  The  physical  characteristics  of  an  AlGaN  layer 
grown  on  top  of  a  GaN  layer  causes  the  production  of  dipoles  in  both  layers.  The 
software  modifications  use  that  concept. 
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Figure  a,  ATLAS-generated  representation  of  AlGaN/GaN  HEMT 

Simulations  were  run  on  a  model  based  on  an  actual  device  tested  at  the  Naval 
Research  Laboratory  (NRL).  The  model  proved  capable  of  matching  the  expected  energy 
band  profile  and  elecli'on  concentration  near  the  heterojunction. 


XIV 


The  model  with  a  matching  Al  mole  fraction  (x  *  0.3)  was  run  through  gate 
biasing  conditions  and  compared  with  measured  results  provided  by  NRL  and  simulated 
previously  by  the  Naval  Postgraduate  School.  The  device  was  simulated  under  matching 
gate  biases  and  drain  voltage  ramping.  The  resultant  IV  curves  displayed  an  above- 
average  correspondence  between  modeled  and  measured  results.  Numerous 
combinations  of  electron  mobility,  velocity  saturation,  and  gale  work  functions  were 
attempted  to  improve  IV  curve  correlation.  The  reported  results  were  derived  from  a 
simulation  run  of  the  reported  modeling  program.  Possible  methods  for  improving  model 
accuracy  are  discussed  and  recommended  for  future  research  and  study. 
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L  INTRODUCTION 


A.  BACKGROUND 

The  advancement  of  radar  and  communication  technology  is  both  vital  and 
necessary  for  the  continued  accomplishment  of  the  Navy's  mission  to:  (1)  protect  friendly 
naval  forces;  (2)  project  national  power;  and  (3)  protect  national  interests.  One  of  the  key 
components  is  the  further  development  of  radio  frequency  (RF)  amplifiers.  Normally,  RF 
waves  are  formed  by  relatively  small  and  physically  fragile  electronics  and,  as  a  result, 
RF  amplifiers  are  used  to  magnify  the  respective  waves  prior  to  transmission.  Many 
systems,  such  as  the  Aegis  Weapons  Systems,  use  travelling-wave  vacuum  tube  (TWT) 
technology  to  perform  the  RF  amplification.  Unfortunately,  TWT's  have  continually 
demonstrated  a  high  level  of  performance  degradation  in  respect  to  service  time  and  a 
high  cost. 

Solid-state  technology,  namely,  semiconductor  amplifiers,  offers  the  most 
promising  alternative  to  TWT  amplifiers.  Incorporating  a  smaller  size  coupled  with 
higher  reliability,  semiconductor  amplifiers  offer  a  greater  ease  of  manufacturing, 
packaging,  and  handling  than  the  TWT's.  Although,  semiconductor  amplifiers  cannot 
currently  match  the  maximum  amplification  of  TWTs,  they  demonstrate  respectable 
power  levels  with  lower  distortion  levels  than  TWT's,  especially  over  extended  periods 
of  time  [1],  [2]. 

Field  Effect  Transistors  (FET's)  and,  specifically.  High  Electron  Mobility 
Transistors  (HEMT's)  are  excellent  RF  power  amplifiers  especially  above  4  GHz  and 
where  low  noise  is  vital.  Due  to  environmental  and  man-made  conditions,  low  noise  is 
essential  for  radar  systems'  optimal  performance.  Although,  HEMT's  are  relatively  new, 
in  relation  to  FET’s,  they  have  demonstrated  great  potential  as  high-speed  power 
amplifiers.  In  general,  FET's  have  displayed  superior  power  gains,  lower  noise,  and 
higher  output  power  capabilities  than  their  historic  competitor,  the  bipolar  junction 
transistor  (BJT)  [3].  However,  the  HEMT,  has  demonstrated  the  highest  current-gain 
cutoff  frequency  of  any  current  semiconductor  device  at  472  GHz  [4]. 
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As  a  compound  semiconductor,  Gallium  Nitride  (GaN)  has  the  potential  to  be  the 
leader  in  RF  amplification  because  of  its  high  frequency  and  high  temperature 
capabilities  and  extremely  high  power  output  [1].  According  to  the  Office  of  Naval 
Research,  GaN  amplifiers  will  take  the  place  of  vacuum  tubes  in  most  Navy  radar 
systems  once  they  are  available.  Because  GaN  can  deliver  up  to  ten  times  more  power  at 
microwave  frequencies  than  current  Silicon  (Si)  or  Gallium  Arsenide  (GaAs) 
semiconductors  used  in  some  radar  systems,  it  is  expected  to  deliver  large  improvements 
in  performance. 

Computer  software  modeling  has  proven  to  be  a  versatile  and  valuable  tool  for 
engineering  design  and  analysis.  Recently,  the  Silvaco  software  applications  have  come 
to  the  forefront  of  engineering  design  and  analysis  of  semiconductor  devices  and 
processes.  In  lieu  of  this,  GaN-based  devices  require  special  attention  and  consideration 
for  computer  modeling  and  simulation  due  to  their  inherent  polarization,  piezoelectric 
(PZ)  and  spontaneous,  properties.  Considerable  research  time  and  funding  has  been 
devoted  toward  these  and  other  nitride-based  devices,  and  has  included  the  establishment 
of  select  constants  for  their  PZ  properties  and  related  effects.  The  polarization  effects  of 
GaN-based  semiconductor  devices  have  not  been  accurately  modeled  in  a  Technology 
Computer  Aided  Design  (TCAD)  program. 

The  GaN's  PZ  effect  is  its  physical  property  to  exhibit  polarization  when  strained. 
This  polarization  results  in  lowered  resistivity  due  to  the  presence  of  extra  charge  carriers 
thus  greater  current  flow  and  power.  Adding  mechanical  strain  to  a  device  might  further 
enhance  PZ  effect  [1]. 

B.  RELATED  WORK 

The  Naval  Research  Laboratory  (NRL)  and  ONR  continue  to  fund  research 
through  educational  institutions  and  industry  on  GaN  HEMT  power  amplifiers.  NRL  is 
modeling  devices  of  their  own  construction  [5].  NRL's  models  do  not  automatically 
account  for  the  piezoelectric  effect.  The  intent  of  this  work  is  to  aid  their  efforts  and  take 
a  further  step  in  accurately  replicating  the  AlGaN/GaN  HEMT. 


2 


The  University  of  California  at  Santa  Barbara  is  computer  modeling  AlGaN/GaN 
HEMT’s  based  upon  the  field  plate  technique  and  surface  donor-like  traps  and  is  focused 
on  breakdown  voltage  research  [6].  Their  models  do  not  automatically  account  for  the 
piezoelectric  effect,  either.  This  work  can  aid  in  developing  closer  correlation  between 
computer  modeling  and  experimental  data. 

The  Naval  Postgraduate  School  has  modeled  AlGaN/GaN  HEMT's  using  the 
Silvaco  Software.  That  work  focused  upon  modeling  an  existing  HEMT's  using 
Silvaco's  C-INTERPRETER  [1].  This  work  will  aid  in  refining  the  previous  modeling 
by  establishing  a  base  Silvaco  code,  using  a  dipole  technique,  to  compare  against 
experimental  data  and  previous  modeling. 

C.  OBJECTIVES 

This  work  seeks  to  increase  the  Navy's  capability  to  computer  model  GaN 
HEMT's  by  providing  working  Silvaco  Software  models  of  an  existing  GaN  HEMT 
device.  The  models  developed  are  distinguished  in  that  they  incorporate  the  piezoelectric 
effect  using  accepted  constants,  a  dipole-concept  approach,  and  computer  program 
generated  values. 
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II.  TRANSISTOR  FUNDAMENTALS 


A.  FIELD  EFFECT  TRANSISTOR 

FET’s  were  the  first  solid-state  transistors,  or  transfer  resistors.  Since  their 
creation  in  the  1950’s,  FET’s  have  distinguished  themselves  from  other  transistors 
through  their  ease  of  production  and  compact  nature.  Thanks  to  the  computer  industry, 
the  FET  (specifically  the  MOSFET)  is  currently  the  most  produced  manmade  object  in 
the  world  [1][7]. 

As  noted  by  the  name,  FET's  operate  based  upon  the  field  effect  phenomenon. 
The  basic  FET  requires  three  contacts  and  a  semiconductor  material.  The  three  contacts, 
source,  gate,  and  drain,  perform  functions  to  facilitate  the  flow  of  electrons  and  current. 
The  source,  an  ohmic  contact,  is  the  electrons’  entry  point  to  the  transistor  from  the 
external  circuit.  The  drain,  an  ohmic  contact,  is  the  electrons'  exit  point  from  the 
transistor  to  the  external  circuit.  The  gate,  a  rectifying  contact,  is  the  “electron  flow" 
controller.  The  gate  controls  electron  flow  by  modulating  the  conductance/resistance  of 
the  semiconductor  material  by  applying  a  voltage  to  the  semiconductor  and  thereby 
modulates  the  current  flow  through  the  FET  (Figure  1).  Electric  current  entering/leaving 
the  FET  is  allowed  to  flow  freely,  restricted,  or  not  at  all  depending  on  the  voltage 
applied  at  the  gate. 

The  gate  to  drain  voltage  difference  will  cause  a  depletion  region  (or  an  inversion 
layer)  to  grow  in  the  vicinity  of  the  gate.  When  the  gate  is  oppositely  biased  from  the 
source-drain  voltage  difference,  the  depletion  region  will  grow  until  it  “pinches  off’ 
current  flow.  As  a  result  of  pinch  off,  the  amount  of  current  flowing  from  source  to  drain 
is  limited  to  its  existing  level.  The  device  will  shut  off  when  no  source-drain  voltage 
difference  exists  and  the  gate  is  biased  enough  to  create  a  complete  depletion  region.  The 
transistor  works  because  the  gate's  voltage  controls  the  size  of  the  depletion  region 
present  in  the  semiconductor. 
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Figure  1.  Basic  FET  Operation. 


B.  HIGH  ELECTRON  MOBILITY  TRANSISTOR 

The  HEMT  has  gone  by  several  different  names:  modulation-doped  field-effect 
transistor  (MODFET);  heterostructure  field-effect  transistor  (HFET);  and  selectively 

doped  heterojunction  transistor  (SDHT).  Regardless  of  the  name  used,  these  devices 
have  unique  performance  characteristics  that  make  them  very  attractive  for  military  and 
civilian  wireless  communications  and  military  radar  applications.  These  characteristics 
are  directly  related  to  the  type  of  carrier  channel  developed  at  the  transistors' 
heterojunction.  The  carrier  channel,  a  quantum  well  with  a  two-dimensional  electron  gas 
(2DEG),  allows  electrons  to  flow  in  only  two  dimensions  vice  three. 

The  2DEG  is  formed  at  the  heterojunction  as  a  result  of  the  conduction  band 

offset  (AEc)  between  the  two  materials.  At  the  interface,  the  material  with  the  narrower 

bandgap  will  have  its  conduction  band  pulled  down  to  the  fermi  level  which  will  make 

the  junction  highly  conductive.  This  material's  conduction  band  rises  with  depth  relative 

to  the  heterojunction  and  guides  electrons  toward  the  material  interface.  On  the  other 

side  of  the  junction,  the  conduction  band  is  much  higher  and  prevents  electrons  from 
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traveling  past  the  material  interface  (Figures  2  &  3).  As  a  result,  a  quantum  well  is 
produced. 


Figure  3.  Three-dimensional  energy  (conduction)  band  diagram  for  unbiased  HEMT. 


In  the  quantum  well,  no  electron  can  be  placed  at  energies  lower  than  the  energy 
level  that  corresponds  to  the  half-wavelength  of  the  electron  [8].  This  energy  level  is 
called  the  first  energy  subband  and  is  below  the  fermi  level  (Ef).  The  second  energy 
subband  corresponds  to  the  electron’s  full-wavelength  and  is  above  Ef.  Due  to  the  slope 
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and  height  of  the  conduction  bands  relative  to  Ef  and  the  first  energy  subband,  electrons 
can  not  move  freely  or  randomly  scatter  in  the  direction  of  the  quantum  well's  cross 
section  (x-axis).  In  Figure  3,  observation  of  the  first  energy  subband  would  suggest 
approximately  zero  electron  movement  in  the  x-direction  but  considerable  movement  in 
the  y  and  z  directions  (two-dimensions).  Relating  the  motion  of  molecules  in  gases  to  the 
motion  of  electrons  in  a  semiconductor  crystal  has  fostered  the  ‘‘electron  gas"  concept. 
The  combination  of  two-dimensional  movement  and  the  “electron  gas"  concept  has 
resulted  in  these  electrons  being  referred  to  as  the  2 DEG. 

The  advantages  of  the  2 DEG  are  that  there  is  reduced  scattering  in  two- 
dimensional  motion  vice  three,  and  this  causes  reduced  noise  levels  and  increased 
electron  mobility.  Both  properties  are  very  desirable  in  high  frequency  and  low  noise 
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applications  such  as,  communications  and  radar.  For  some  current  applications,  HEMTs 
have  been  constructed  of  non-Silicon  semiconductors  such  as,  GaAs  and  GaN,  due  to  the 
specific  material  properties  of  these  compound  semiconductors  (Table  1).  Additionally, 
for  AlGaN/GaN,  2DEG  mobility  was  found  to  be  1.9x10'*  cm' /Vs  @  lOK  [9]. 

C.  GALLIUM  NITRIDE 

Group  Ill-nitride  materials  arc  ideal  for  high  power  and  high  temperature  devices 
due  to  their  large  energy  bandgap,  high  breakdown  voltage,  high  peak  electron  velocity 
and  high  electron  sheet  density  in  channels  when  used  in  a  hetero structure  [1].  Nitride 
materials  benefit  from  having  large  and  direct  bandgaps  which  makes  them  very  capable 
and  well  suited  to  handle  environmentally  hostile  conditions  such  as,  those  at  high 
altitudes  or  in  space,  and  very  useful  in  optical  applications.  Although  other  materials 
such  as,  Silicon  Carbide  (SiC),  have  some  similar  properties,  nitride  materials  can  form 
better  hetero  structures  and  ohmic  contacts,  and  appear  to  be  more  promising.  Of  the 
nitride  materials,  the  GaN  based  devices  have  demonstrated  the  most  functionality  and 
garnered  the  most  interest.  Figure  4  is  a  depiction  of  a  Wurtzite  GaN  molecule. 


Figure  4.  Depiction  of  Wurtzite  GaN  molecule  with  representative  hexagonal  axes 
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The  GaN  devices’  major  drawbacks  are:  (1)  challenges  in  the  manufacturing 
process  that  have  lead  to  reliability  problems  and  excessive  costs;  and  (2)  electron 
mobility  is  not  significant  for  power  devices  at  high  voltages.  Along  with  other  Nitride 
materials,  GaN  has  been  difficult  to  chemical  etch  because  of  its  characteristically  inert 
nature.  Because  of  their  high  melting  temperatures  and  pressures,  normal  etching 
methods  such  as,  the  Czochralski  method,  have  proven  ineffective  on  Nitride  materials. 
The  commercial  semiconductor  industry,  some  educational  institutions,  and  various 
research  entities  have  performed  significant  research  into  the  growth  of  quality  GaN 
crystals  to  produce  better  performance  and  fewer  defects,  however,  further  work  needs  to 
be  performed. 

In  general,  high  temperatures  degrade  electronic  performance  because  electric- 
field  mobility  and  electric-field  electron  velocities  can  be  adversely  affected  by  rising 
temperature.  For  the  production  of  GaN  devices,  like  all  high  power  devices,  the 
problem  of  thermal  dissipation  had  to  be  addressed.  This  has  made  the  search  and 
selection  of  a  suitable  substrate  of  paramount  importance. 

In  the  past.  Sapphire  has  been  a  successful  substrate  for  GaN  devices.  Sapphire 
was  the  first  successful  substrate  used  with  GaN  by  Shui  Nakamura  (who  made  a  blue 
LED)  at  Nichia  Chemical  in  Japan  [10].  Silicon  Carbide,  however,  can  conduct  heat  at 
seven  times  the  rate  of  Sapphire  [11].  Thus,  Silicon  Carbide  appears  to  be  a  better  GaN 
substrate  material  than  Sapphire.  Individually,  SiC  is  a  useful  wide  bandgap 
semiconductor,  as  evidenced  by  its  material  properties  in  Table  1,  and  is  much  easier  to 
fabricate  than  Sapphire.  In  addition  to  boasting  a  much  better  thermal  conductivity  and  a 
greater  ease  of  fabrication,  SiC  has  a  lower  difference  between  lattice  constants  (3% 
difference  for  GaN/SiC  compared  to  13%  for  GaN/Sapphire)  which  reduces  the 
molecular  stress  between  semiconductor  and  substrate.  Also,  as  a  result  of  lower 
dislocation  density,  SiC  produces  a  larger  2 DEG  electron  concentration  (--1.4xlO*Vcm^ 
vice  lxlO*Vcm^)  [12].  Despite  SiC's  advantages.  Sapphire  still  continues  to  be  a  popular 
substrate  because  of  its  significantly  lower  cost.  Another  low  cost  alternative  that  has 
been  researched  is  Si.  Si  substrates  not  only  offer  a  low  cost  alternative  with  roughly  one 
half  the  thermal  conductivity  of  SiC  but  also  allow  the  potential  integration  of  power 
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electronics  on  an  advanced  and  mature  Si  technology  [13].  Other  substrate  candidates 
that  have  been  investigated  include  Lithium  Gallate  [14]  and  Neodymium  Gallate  [15]. 


D.  POLARIZATION  EFFECTS 

Certain  molecular  crystals  become  polarized  when  they  are  mechanically  stressed 
(Figure  5)  [16].  When  the  atoms  are  pushed  from  their  natural  positions,  their  charges  do 
not  balance  appropriately.  As  a  result,  surface  charges  appear  and  a  voltage  difference 
between  the  two  crystal  faces  is  produced.  These  same  crystals  will  also  exhibit 
mechanical  strain  or  distortion  when  they  experience  an  electric  field.  The  direction  of 
distortion  will  depend  upon  the  direction  of  the  electric  field  or  polarity  of  the  applied 
voltage.  These  characteristics  describe  the  PZ  effect  inherent  to  these  crystals. 


Only  crystals,  such  as  group  Ill-nitrides,  without  a  center  of  symmetry,  called 
noncentrosymmetric,  can  exhibit  PZ  properties.  In  addition  to  the  strain-induced  PZ 
polarization,  calculations  have  shown  that  group  lll-nitrides  have  strong  spontaneous 
polarization  [17]  [18].  In  GaN,  it  can  be  seen  that  the  direction  of  polarization  is  along  the 
c-axis  from  N  to  Ga  (Figure  5).  The  direction  of  a  device's  polarization,  Ga-  or  N-faced, 
is  very  important  because  it  will  determine  where  the  2DEG  will  form  in  a  HEMT. 

The  face  type  of  a  GaN-based  device  determines  the  location  of  the  2DEG  but  the 
percentage  of  aluminum  (Al)  determines  the  charge  density  in  A1  GaN /GaN  HEMT's, 
Studies  have  shown  a  direct  relationship  between  the  concentration  of  Al  in  the  AlGaN 
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layer  to  sheet  charge  density  at  the  interface  [19]  because  AlGaN  has  a  different  lattice 
constant  than  GaN  and  in-plane  biaxial  stress  is  created  when  AlGaN  is  grown  on  GaN 
(Figure  6)  [20].  This  is  the  manifestation  of  the  piezoelectric  effect  in  an  AlGaN/GaN 
HEMT. 

GaN  and  its  compounds  have  high  piezoelectric  properties  that  result  in 
measurable  electric  fields  in  the  mega-Volts  per  centimeter  range.  These  high  electric 
fields  effect  conductivity,  allowing  for  electron  transport  without  using  dopants  [1]. 


Figure  6.  Depiction  of  AlxGai.x-faced  material  grown  onto  Ga-faced  material. 


(Negative  charge  is  vertically  up,  positive  charge  is  vertically  down) 


Much  of  the  past  research  into  the  polarization,  PZ  and  spontaneous,  properties 
between  AlGaN  and  GaN  layers  has  focused  on  the  internal  molecular  strain  caused  by 
their  lattice  mismatch  (Figure  7).  Of  course,  the  amount  of  external  stress  or  strain  is 
another  area  of  research  and  factor  to  consider  in  the  accurate  evaluation  of  AlGaN/GaN 
semiconductor  performance. 
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Figure  7.  Bandgaps  and  lattice  constants  for  select  semiconductor  compounds. 


Additionally,  an  applied  electric  field  might  influence  the  PZ  effect.  In  simplistic 
terms,  piezoelectric  crystals  are  electromechanical  transducers  that  change  strain  into 
electric  potential  and  vice  versa.  Theoretically,  the  application  of  an  external  electric 
field  will  cause  some  mechanical/physical  defonnation  that  could  further  affect  the  total 
PZ  effect.  Another  potentially  important  parameter  is  thermal  expansion  strain.  Further 
research  into  developing  reliable  relationships  between  layer  microstructure  and  thermal 
strain  is  needed  to  potentially  include  this  parameter  in  future  modeling. 
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Ill,  DEVICE  MODELING  AND  SIMULATION 


A.  SILVACO 

Silvaco  International  produced  the  device  modeling  and  simulation  software 
utilized  in  this  work.  Silvaco’s  ATLAS™  is  a  versatile  and  modular  program  designed 
for  one,  two,  and  three-dimensional  device  simulation.  BLAZE™  and  GIGA™, 
ATLAS™  sub-modules  (Figiue  8),  perfonn  specialized  functions  required  for  advanced 
materials,  heterojunctions,  and  temperature- dependent  conditions.  To  control,  modify, 
and  display  the  modeling  and  simulation,  the  Virtual  Wafer  Fabrication  (VFW) 
Interactive  Tools,  namely  DECKBUILD™  and  TONYPLOT™,  were  utilized  (Figure  9). 
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Figure  8.  Representation  of  ATLAS’  modular  structure. 
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Unlike  some  other  modeling  software,  Silvaco  uses  physics-based  simulation 
rather  than  empirical  modeling.  In  truth,  empirical  modeling  produces  reliable  formulas 
that  will  match  existing  data  but  physics-based  simulation  predicts  device  performance 
based  upon  physical  stmctiue  and  bias  conditions.  To  perform  the  modeling,  the  Silvaco 
software  graphically  represents  a  device  on  a  two-dimensional  grid  with  designated 
electronic  meshing  parameters.  At  every  mesh  intersection,  the  program  simulates  carrier 
transport  by  means  of  differential  equations  derived  from  Maxwell’s  laws.  To  achieve 
accuracy,  the  program  incorporates  the  appropriate  physics  via  numerical  procedures. 

To  accurately  model  the  II I- V  semiconductors,  ATLAS  must  employ  the  BLAZE 
program  extension  to  modify  calculations  that  involve  energy  bands  at  heterojunctions. 
The  heterojunctions  require  changes  in  calculating  current  densities,  thermionic 
emissions,  velocity  saturation,  and  recombination-generation. 

ATLAS  attempts  to  find  solutions  to  carrier  pai'ameters  such  as  current  through 
electrodes,  carrier  concentrations,  and  electric  fields  throughout  the  device.  ATLAS  sets 
up  the  equations  with  an  initial  guess  for  parameter  values  then  iterates  through 
parameters  to  resolve  discrepancies.  ATLAS  will  alternatively  use  a  decoupled 

(Gummel)  approach  or  a  coupled  (Newton)  approach  to  achieve  an  acceptable 
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correspondence  of  values.  When  convergence  on  acceptable  values  does  not  occur,  the 
program  automatically  reduces  the  iteration  step  size.  ATLAS  generates  the  initial  guess 
for  parameter  values  by  solving  a  zero-bias  condition  based  on  doping  profiles  in  the 
device. 

B.  PRIOR  SILVACO  MODELING  EFFORTS 

Other  efforts  have  been  performed  to  model  AlGaN/GaN  HEMT's.  Three  such 
efforts  have  been  conducted  by  Naiqian  Zhang  [21]  and  Shreepad  Karmalkar  [7]  of  the 
University  of  California  at  Santa  Barbara  (UCSB)  and  Karl  Eimers  [1]  of  the  Naval 
Postgraduate  School  (NFS). 

1.  Surface  Donor-Like  States  Approach 

In  [21]  and  [7],  UCSB  follows  Ibbetson’s  theory  [22]  of  surface  donor-like  states 
being  the  source  of  the  electrons  for  the  HEMT's  2DEG  and,  as  a  result  of  the  strong 
internal  electric  field,  the  electrons  are  driven  into  the  quantum  well.  Using  this  concept 
and  the  simplified  approach  to  the  HEMT's  charge  distribution  in  equation  (1), 

"s  ="/-  =  (1) 

where,  is  the  channel  electron  concentration;  n^is  the  electron  sheet  concentration;  np+ 
and  ftp,  are  polarization  dipole  charges  at  opposite  faces  of  the  AlGaN  layer;  n,  is  the 
insulator-donor  layer  interface  charge;  and  is  the  ionized  unintentional  doping  charge, 
both  [21]  and  [7]  attained  results  comparable  with  experimental  data  for  unintentionally- 
doped  and  field  plate  HEMT's,  respectively.  This  simplified  approach  does  not  consider 
the  spontaneous  and  piezoelectric-based  charges  separately. 

2.  C-Interpreter  Approach 

In  [1],  Eimers  uses  Ambacher’s  [23]  equations  (Table  2)  in  ATLAS'  C- 
INTERPRETER  to  model  the  carriers'  characteristics  in  the  AlxGai.xN/GaN  HEMT  (A1 
mole  fraction  x  =  0.3). 
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EQUATIONS 


Sheet  electron  concentration 

Dielectric  constant 

£(Ar)  =  -0.3Ar  +  10.4 

Work  fiinction 

=(l.3x  +  0.84)eF 

Conduction  band  offset 

A£^=0.7KW-£,,(0)J 

Bandgap 

{x)  =  (a:6.  13  +  (1  -  x)3.42  -  x(\  -  x)\  .0>F 

PZ  effect  polarization 

Lattice  constant 

«(,  (x)  =  {-0.077X  +  3. 1 89)1 0  m 

Elastic  constant 

C,3(Ar)  =  (5x  +  103)C?Pa 

Elastic  constant 

C^^(x)^{-32x  +  405)GPa 

Piezoelectric  coefficients 

(X)  =  (e„  (AIN)  -  e.  (GaN))*x  +  (GaN) 

Spontaneous  Polarization 

P,,(x)  =  i-0.052x  -  0.029)-^ 

m 

Degree  of  Relaxation 

0  0<x<OJ8 

3.5Ar-1.33  0.38  <jc<  0.67 

1  0.67  <  a:  <1 

Sheet  charge 

o'(^)  =  {Pp£  (fyollom)  +  Py,(hoilom)}~  +  ^sp 

Fermi  level 

(x)  =  -0. 1 02967  +  {2.191 7).*  -  0}x^ 

Table  2.  Equations  utilized  in  [1]. 
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Eimers  concentrated  on  utilizing  these  equations  in  ATLAS’  C-INTERPRETER  in 
conjunction  with  a  DEVEDIT-generated  HEMT  model  to  simulate  the  polarization 
effects  and  electron  sheet  concentration  within  the  transistor.  Despite  DEVEDlT’s 
inability  to  recognize  AlGaN  as  a  material,  Eimers  was  able  to  modify  the  program’s 
material  properties  for  AIN  to  approximate  AlGaN.  Figure  10  is  a  representation  of  the 
resultant  IV  curves  for  OV,  -IV,  -3V,  and  -5V  gate  bias  with  drain  voltage  being  ramped 
from  OV  to  20V.  Figure  11  displays  the  experimentally  derived  IV  curves  for  the  doped 
AlGaN/GaN  HEMT. 


Figure  10.  Simulated  IV  curves  for  doped  AlGaN/GaN  HEMT.  (From:  [1]) 
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Figure  1  L  Experimentally  derived  IV  curves  for  doped  AlGaN/GaN  HEMT.  (From:  [1]) 


The  modeled  IV  ciures  had  noticeable  differences.  Despite  the  differences,  the 
modeled  and  experimental  curves  were  roughly  similar.  This  demonstrated  the  potential 
for  using  the  Silvaco  software  installed  on  a  UNIX  system,  specifically  the  C- 
INTERPRETER,  to  model  doped  HEMT’s. 

C.  THIS  EFFORT 

This  work  focuses  on  an  AlosGao  7N/GaN  HEMT  that  was  constructed  and 
physically  tested  by  S.  C.  Binari  (et  al.)  at  the  Naval  Research  Laboratory  (NRL)  in 
Washington,  D.C.  and  is  the  second  approach  by  NFS  to  model  this  HEMT.  Figure  12  is 
a  representation  of  a  HEMT  tested  in  [24]  and  Figure  13  is  the  ATLAS-generated 
representation  of  the  HEMT. 

Noted,  the  model  appears  different  from  the  tested  HEMT.  First,  the  source  and 
drain  contacts  were  placed  on  the  sides  of  the  AlGaN  and  GaN  layers  to  allow  proper 
current  flow  from  contacts  to  the  device,  The  change  in  the  contact  placement  is  required 
because  the  software  treats  the  AlGaN  layer  as  an  effective  insulator  preventing  any 
current  flow.  In  addition,  the  implant  damage  zones,  which  provide  isolation  in  the  real 
device,  were  not  included. 
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Figure  12.  Representation  of  NRL’s  doped  AlGaN/GaN  HEMT. 


Figure  13.  AT  LAS -gen  crated  representation  of  NRL*s  doped  AlGaN/GaN  HEMT 


21 


For  modeling  purposes,  this  work  proceeded  using  an  approach  aimed  at 
representing  the  HEMT  layers,  AlGaN  and  GaN,  as  dipoles.  With  this  dipole  approach, 
we  assumed  Ga-faced  and  AlGa- faced  layers,  respectively.  To  simulate  the  HEMT’s 
characteristic  2DEG  and  polarization,  the  AlGaN  layer  is  n-doped  5e+18  from  the  top  of 
the  layer  to  0.005  microns  then  p-doped  5e+18  from  0.019  to  0.024  microns,  and  the  GaN 
layer  is  n-doped  2e+19  from  0.030  to  0.035  microns  then  p-doped  2e-i-18  from  3.0  to 
3.024  microns  (Figure  14,  Appendix  A).  As  can  be  observed,  the  top  layer  of  electron 
concentration  width  in  the  AlGaN  layer  is  approximately  as  entered  but  both  the  bottom 
layer  of  AlGaN  hole  concentration  width  and  top  layer  of  GaN  electron  concentration 
width  are  decreased  due  to  recombination.  Deckbuild’s  default  recombination  lifetimes, 
1x10*’  sec,  were  utilized  for  this  effort.  Taking  a  cutline  across  this  section  (vertically 


Figure  14.  TONYPLOT’s  representation  of  doping  profile. 


down  from  the  the  below  electron  concentration  profile  is  produced  (Figure  15). 
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Figure  15.  Vertical  cutline  through  AlGaN  layer  and  top  of  GaN  layer. 


Taking  peak  data  from  Figure  15  and  using  an  algebraic  summation  similar  to  equation 
(1),  the  result  is: 

(5xl0‘*  donors/cm^)-  (5x10*^  acceptors/cm")  +  {IxlO'^  donors/cm^)  =  IxlO’^  donors/cm* 


This  very  simple  computation  does  not  account  for  the  bottom  of  the  GaN  layer  because 
its  interaction  with  the  AlGaN  layer  and  top  of  the  GaN  layer  can  be  assumed  to  be 
negligible  (3  fim  vice  0.025  [im). 
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IV.  RESULTS 


A.  2-DEG 

Making  a  ciitline  from  the  gate  at  zero  bias  to  the  AIN  layer.  Figure  16  displays 
the  desired  spike  of  electrons  on  the  GaN  side  of  the  AlGaN/GaN  interface.  The  dark 
line  is  the  interface  with  AlGaN  on  the  left  and  GaN  on  the  right.  The  spike  in  electron 
concentration  (“-0.8e+l9/cm^)  is  representative  of  the  2-DEG  resulting  from  polarization. 
Figure  17  is  the  conduction  band  energy  diagram  for  zero  gate  bias.  The  “dip"  or 
quantum  well  drops  to  a  -0.05  eV  energy  level  to  facilitate  electron  transport.  The 
HEMT’s  characteristic  quantum  well  reduces  electron  scattering  and  increases  electron 
mobility. 


Figure  16.  Electron  concentration  at  the  AlGaN/GaN  interface  with  zero  gate  bias. 
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B.  ELECTRON  FLOW 

Figures  18-21  are  samples  of  TONYPLOT’s  representations  of  the  electron 
concentrations  in  the  AlGaN/GaN  HEMT  under  OV  gate  bias  and  drain  voltages  OV  and 
20V.  Graphics  for  gate  biases  -IV,  -2V,  -3V,  -4V,  and  -5V  are  in  Appendix  B.  These 
displays  will  indicate  the  flow  of  electrons  (and  holes)  through  the  HEMT.  Figures  18 
and  19  are  the  zero  gate  bias  (Vg)  and  OV  drain  voltage  (Vd)  condition  displays.  As  can 
be  see,  the  electron  concentration  is  highest  near  the  AlGaN/GaN  interface.  In  Figure  19, 
it  is  clearly  shown  how  the  ohmic  gate  negates  the  electron  concentration.  Figures  20  and 
21  are  the  Vg  =  OV,  V^  =  20V  condition  graphics.  A  decreasing  electron  concentration 
can  be  observed  below  the  right  side  of  the  gate.  This  is  showing  the  development  of  the 
transistor's  depletion  region.  Figures  22-26  are  samples  of  TONYPLOT's 
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Figure  18.  HEMT  under  zero  gate  bias  and  OV  drain  voltage. 


Air  us 

Daca  fcoa  HP5Ga}(Ei3BitO  Q.str 


'  '  '  I  '  '  '  J  '  •  i  •  •  •  I  '  '  '  I  '  '  I  '  I  I 

2.4  3X  2.0  i  3A  2.4  3.0  3,2  4 

Hierara 


Figure  19.  Close-up  view  of  HEMT  under  zero  gate  bias  and  OV  drain  voltage. 


27 


ATLAS 

Data  from  WSGaNHeratO  l.sCc 


0  1  3  3  4  S  6 

Mictorts 


Figure  20.  HEMT  under  Vg  =  OV,  Vj  =  20V  condition.  Depletion  region  developing. 
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representations  of  the  current  density  in  the  AlGaN/GaN  interface  under  OV  gate  bias  and 
drain  voltages  OV,  5V,  lOV,  15V,  and  20V.  As  can  be  observed  in  Figure  22,  the  region 
has  a  positive  current  density  (--3790  A/cm").  However,  after  drain  voltage  increases 
(Figure  23),  a  negative  current  density  develops  below  the  interface.  This  negative 
current  density  represents  the  2DEG  and  negative  current  flow.  With  increasing  drain 
voltage  (Figures  23-26),  a  decrease  and  disniption  of  the  negative  ciu'rent  flow  is  shown, 
The  disruption  or  “pinch  off'  is  caused  by  the  development  of  the  depiction  region. 
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Figure  23,  Close-up  of  HEMT  under  Vg  =  OV,  Vd  =  5V  condition. 
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Figure  25.  Close-up  ofHEMT  under  Vg  =  OV,  Vj  =  15V  condition. 


-0.2  - 


-0.15  — 

-□.1  — 

-0.05  — 


0 

0.05 

u.l 


Microns 


Figure  26.  Close-up  of  HEMT  under  Vg  =  OV,  Vd  =  20V  condition. 
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c. 


IV  CURVES 


The  Alo.3Gao.7N/GaN  HEMT  model  was  tested  and  attained  the  below  IV  curves 
for  Vg  =  OV,  -IV,  -2V,  -3V,  -4V,  and  -5V  while  ramping  Vd  from  0-20V  (Figure  27). 
This  work's  results  were  compared  with  NRL's  measured  results  [24] (Figure  29)  and 
Eimers'  simulated  results  [1]  (Figure  30)  for  Vg  =  OV,  -IV,  -3V,  and  -5V  (Figure  28). 


ATLAS  OVERLAY 
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'igure  27.  Simulation-generated  IV  curves  for  OV,  -IV,  -2V,  -3V,  -4V,  and  -5V. 


In  Figures  28  and  29,  the  *‘knee'’  voltages  are  corresponding  (Vd  ^  3V)  for  the  OV  and  - 
IV  ciuves.  In  addition,  the  slope  of  the  OV,  -IV,  -3V  curves  (:s0.2  mA/V)  show  strong 
correlation.  The  dissimilarities  noted  are  the  “shortfalls"  on  the  curve  peaks  and  the 
current  collapse.  The  peak  “shortfalls"  appear  to  be  a  result  of  the  depletion  region’s 
growth  in  the  simulation.  Multiple  modifications  to  the  gate's  work  function  did  not 
resolve  this  occurrence.  The  current  collapse  in  the  experimental  curves  was  theorized  to 
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be  a  result  of  GaN  defect  traps  [24].  Binari  considered  the  graphically-displayed  current 
collapse  to  be  consistent  with  the  observed  correlation  between  the  amount  of  current 
collapse  and  resistivity  of  the  buffer  layer.  Defect  traps  were  not  simulated  in  tliis  work 
thus  the  program  was  unable  to  compensate  for  these  potential  effects.  Figure  31  is  an 
overlap  of  NRL’s  and  this  work's  IV  curves  for  Vg  =  OV,  -IV,  -3V,  and  -5  V. 
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Figure  30.  Eimers’  simulated  IV  curves  for  OV,  -IV,  -3V,  and  — 5V  gate  bias.  (From:  [1]) 
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Figure  31.  Overlap  ofNRL’s  and  this  work's  IV  curves  for  Vg  =  OV,  -IV,  -3V,  and-5V. 


Appendix  C  contains  ti'ans conductance  and  subthreshold  plots  of  the  Alo,3Gao,7N/GaN 
HEMT.  As  can  be  observed  in  Figure  Cl,  threshold  voltage  is  approximately  —4,5  volts. 
Experimental  curves  of  these  were  not  available  thus  no  comparison  was  performed. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

The  polarization,  piezoelectric  and  spontaneous,  effects  are  significant  in 
AlGaN/GaN  devices  and  can  be  modeled  with  some  degree  of  accuracy  utilizing  the 
Silvaco  software.  In  this  work,  the  desired  density  of  carriers  has  been  demonstrated  to 
concur  with  established  theory  by  performing  the  modeling  using  a  dipole  approach. 
Also,  the  current  versus  voltage  performance  (IV  curves)  of  the  modeled  device 
approximates  a  measured  device  to  a  degree  of  accuracy. 

Nevertheless,  the  model  requires  greater  refinement  and  treatment  to  more  closely 
match  actual  device  performance.  Some  potential  areas  for  cause  and  correction  of  the 
inaccuracies  have  been  previously  outlined.  In  addition  to  these,  other  potential  methods 
to  resolve  the  IV  curve  discrepancies  are  currently  available  in  the  Silvaco  software. 
Some  of  methods  include  modeling  the  thermal,  interface,  and  quantum  effects  related  to 
the  device  and  its  environment.  The  Silvaco  software  has  a  sub- program,  GIGA™, 
which  extends  ATLAS^^  to  account  for  heat  transfer  at  the  lattice-level  by  implementing 
Wachutka's  thermodynamic  model.  Also,  ATLAS™  has  an  INTERFACE  statement  that 
allows  one  to  define  the  interface  charge  density.  This  function  might  allow  for  a  simple 
method  for  defining  the  2-DEG  but  it  could  modify  the  surface  recombination  velocity 
and  thermionic  emissions,  which  might  be  undesirable.  The  quantum  effects  can  be 
addressed  in  Silvaco  by  solving  Schrodinger's  equation,  which  will  modify  the  normally 
calculated  density  of  states  and  carrier  concentrations. 

This  work  reported  mainly  the  results  derived  from  a  ATLAS  program  that  had 
the  piezoelectric  and  spontaneous  polarizations  entered  as  a  combination  (Appendix  A). 
However,  the  same  results  were  attained  using  a  program  which  algebraically  separated 
the  polarizations  (Appendix  D). 
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B.  RECOMMENDATIONS 

Strenuous  and  concentrated  efforts  should  be  placed  on  modeling  a  closer  match 
to  real  device  performance.  Revisit  the  C-INTERPRETER  approach  to  modeling  using 
some  lessons  learned  from  this  work  because,  once  it  is  refined,  the  C-lNTERPRETER 
method  should  provide  greater  flexibility  for  various  HEMT's  and  conditions. 

Work  closer  with  Silvaco  to  incorporate  the  piezoelectric  and  spontaneous 
polarization  effects  into  their  modeling  software.  This  would  aid  the  Navy  in  researching 
and  developing  newer  GaN  devices.  Additionally,  investigation  into  the  possibility  of 
modeling  the  dynamic  effects  between  the  external  stress,  applied  voltage  bias,  and 
polarization  effects  should  occur. 

For  Naval  Postgraduate  School  (NPS)  Electrical  Engineering  students,  an 
introduction  to  Silvaco 's  modeling  software  at  the  NPS  undergraduate-course  level.  This 
would  provide  students  and  staff  with  additional  corporate  knowledge  on  a  software 
package  that  is  currently  being  employed  by  sponsors,  such  as  the  Navy  Research 
Laboratory,  the  Space  and  Warfare  Command,  and  other  institutions  involved  with 
semiconductor  research. 
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APPENDIX  A.  DECKBUILD  SIMULATION  INPUT  DECK 


go  atlas 

Title  AlGaN/GaN  HEMT  with  AIN  Nucleation  Layer 

#  based  upon  NRL  experimental  HEMT  (0.3  Al  mole  fraction) 

#  Vg  =  0V 

# 

# 

#  SECTION  1 :  Meshing  Parameters 

# 

mesh 


x.mesh  loc=0.0 

spac=1.0 

x.mesh  loc=1.0 

spac=0.1 

x.mesh  loc=3.0 

spac=0.05 

x.mesh  loc=5.0 

spac=0.05 

x.mesh  loc=6.0 

spac=1.5 

y.mesh  loc=0.0 

spac=0.005 

y.mesh  loc=0.025  spac=0.005 

y.mesh  loc=1.0 

spac=0.05 

y.mesh  loc=3. 025  spac=2.0 

y.mesh  loc=3.055  spac=2.0 

y.mesh  loc=8.0 

spac=3.0 

# 

# 

#  SECTION  2:  Structure  Specification 

# 

region  num=I  material=AlGaN  y.min=0.0  y.max=0.025  x.composition=0.3 
region  num=2  material=GaN  y.min=0.025  y.max=3.025 
region  num=3  materia  1= AIN  y.min=3.025  y.max=3.045 
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region  num=4  materia  1= Sapphire  y.min=3.045  y.max=8 

# 

elec  num=  1  name=source  x.min=0.0  x.max=0.0  y.min=0.0  y.max=  1 .0 
elec  num=2  name=gate  x.min=2.5  x.max=3.7  y.min=-0.5  y.max=0.0 
elec  num=3  name=drain  x.min=6.0  x.max=6.0  y.min=0.0  y.inax=1.0 
# 

doping  uniform  y.min=0.0  y.max=0.005  n.type  conc=5.el8 
doping  uniform  y.min=0.0l9  y.max=0.024  p-type  conc=5.el8 
doping  uniform  y.min=0.03  y.max=0.035  n.type  conc=2.el9 
doping  uniform  y.min=3.000  y.max=3.024  p-type  conc=2.el8 

# 

# 

# 

#  SECTION  3:  Material  Models 

# 

material  mated al=AlGaN  mun=600  mup=10  affinity=3.82  nc300=2.07el8 
material  material=AlGaN  nv300=I.16el9  eg300=3.96  align=0.8 

p  ermitti  vity  =  1 0 . 3  2 

# 

material  material=GaN  mun=900  mup=  1 0  nc300=  L07e  1 8 
material  material=GaN  nv300=  1 . 1 6e  1 9  eg300=3 .4  align=0. 8 
material  material=GaN  arichn=24  arichp=96  edb=0.025  eab=0. 160 
material  material=GaN  permittivity=9.5  vsatn=2e7 
# 

model  material=AlGaN  conmob  fldmob  srh  print 
model  material=GaN  conmob  fldmob  srh  print 
# 

impact  bn=3.4e7  an=2.9e8 
impact  bn=3.4e7  an=2.9e8 
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# 

contact  name=gate  workfun=5.0 

# 

#  SECTION  4:  Id-Vd  calculation 

# 

method  gummel  newton  itlim=20  trap  maxtrap=10  vsatmod.inc=0.01  carriers=2 
output  con.band  val.band  band.param  e. mobility  h. mobility 
solve  vgate=0 
# 

save  outf=NPSGaNHemtO_0.str 
tony  plot  NPSGaNHemt0_0.slr 
# 

log  outf=NPSGaNHemt0_0.1og  master 

# 

method  newton  trap  itlim=3  5  maxlrap=6  carriers=2 
solve  vdrain=0.0  vstep=0.1  name=drain  vfmal=l  vsource=0.0  vgate=0.0 
solve  vdrain=  1.0  vstep=0.5  name=drain  vfmal=5  vsource=0.0  vgate=0.0 
solve  vdrain=5.0  vstep=0.75  name=drain  vfinal=20  vsource=0.0  vgate=0.0 
# 

save  outf=NPSGaNHemtO_Lstr 
tonyplot  NPSGaNHemt0_0.log 
# 

quit 
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APPENDIX  B,  ELECTRON  CONCENTRATION  GRAPHICS 


1.  0  VOLT  GATE  BIAS 


Figure  B 1 .  Vg  =0  volts,  Vd  =  0  volts 


Figure  B2.  Close-up  of  Vg  =0  volts,  Vd  =  0  volts 
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Figure  B3.  Vg  =0  volts  Vd  =  20  volts 


Figure  B4.  Vg  =0  volts  Vd  =  20  volts 
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2 


1  VOLT  GATE  BIAS 


Figure  B5.  Close-up  Vg  =  volt,  Vd  =  20  volts 
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3 


2  VOLT  GATE  BIAS 


Figure  B6.  Vg  =  -2  volts,  =  20  volts 


Figure  B7.  Close-up  Vg  =  -2  volts,  Vd  =  20  volts 
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4 


3  VOLT  GATE  BIAS 


Figure  B8.  Vg  =  -3  volts,  Vd  =  0  volts 


Figure  B9.  Vg  =  -3  volts,  Vd  =  20  volts 
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5. 


•^4  VOLT  GATE  BIAS 


Figure  B 1 1 .  Close-up  Vg  =  -4  volts,  Vj  =  20  volts 


51 


6. 


•5  VOLT  GATE  BIAS 


Figure  B 12.  Vg  = -5  volts,  Vd  =  0  volts 
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MIAS 

Data  fcon  NPSCaMHeatS  l.atr 


Kiesasc 

Figure  B 14.  Close-up  Vg  =  -5  volts,  Vj  =  20  volts 
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APPENDIX  C.  TRANSCONDUCTANCE  AND  SUBTHRESHOLD 


Voltaoe  CV) 

Figured.  Transconductance  Plot 


Sate  voltage  (V) 


Figure  C2.  Subthreshold  Plot 
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APPENDIX  D.  DECKBUILD  SIMULATION  INPUT  DECK 


go  atlas 

Title  AlGaN/GaN  HEMT  with  AIN  Nucleation  Layer 

#  based  upon  NRL  experimental  HEMT  (0.3  Al  mole  fraction) 

#  Vg  =  OV 

#  Simulated  separation  of  charges 

# 

# 

#  SECTION  I:  Meshing  Parameters 

# 

mesh 

x.mesh  loc=0.0  spac=1.0 
x.mesh  loc=  1 .0  spac=0. 1 
x.mesh  loc=3.0  spac=0.05 
x.mesh  loc=5.0  spac=0.05 

x. mesh  loc=6.0  spac=1.5 

y. mesh  loc=0.0  spac=0.005 
y.mesh  loc=0.025  spac=0.005 
y.mesh  loc=  1 .0  spac=0.05 
y.mesh  loc=3. 025  spac=2.0 
y.mesh  loc=3.055  spac=2.0 
y.mesh  loc=8.0  spac=3.0 

# 

# 

#  SECTION  2:  Structure  Specification 

# 

reg  i  0  n  num=  I  materi  a  1= AIG  aN  y .  min=0 . 0  y .  max=0 . 02  5  x .  compo  s  i  tion=0 . 3 
region  num=2  materi al=GaN  y.min=0.025  y.max=3.025 


57 


region  num=3  material=AlN  y.niin=3.025  y.max=3.045 
region  nuni=4  materia  1= Sapphire  y.min=3.045  y.max=8 
# 

elec  num=  1  name=source  x.min=0.0  x.max=0.0  y.min=0.0  y.max=  1 .0 
elec  num=2  name=gate  x.min=2.5  x.max=3.7  y.min=-0.5  y.max=0.0 
elec  num=3  name=drain  x.min=6.0  x.max=6.0  y.min=0.0  y.inax=1.0 
# 

#  Separate  Charges 

#  Simulate  Spontaneous  Polarization 

doping  uniform  y.min=0.0  y.max=0.005  n.type  conc=2.5el8 
doping  uniform  y.min=0.0l9  y.max=0.024  p-type  conc=2.5el8 
doping  uniform  y. mi n=0. 03  y.max=0.035  n.type  conc=l.el9 
doping  uniform  y.min=3. 000  y.max= 3. 024  p-type  conc=l.el8 

#  Simulate  Piezoelectric  Polarization 

doping  uniform  y.min=0.0  y.max=0.005  n.type  conc=2.5el8 
doping  uniform  y.min=0.0 19  y.max=0.024  p-type  conc=2.5el8 
doping  uniform  y.min=0. 03  y.max=0.035  n.type  conc=l.el9 
doping  uniform  y.min=3.000  y.max=3.024  p-type  conc=l.el8 

# 

# 

# 

#  SECTION  3:  Material  Models 

# 

material  material=AlGaN  mun=600  mup=10  affinity=3.82  nc300=2.07el8 
material  material=AlGaN  nv300=I.I6el9  eg300=3.96  align=0.8 

p  ermitti  vity  =  1 0 . 3  2 

# 

material  material=GaN  mun=900  mup=  10  nc300=  1 .07e  1 8 
material  material=GaN  nv300=  1 . 1 6e  1 9  eg3  00=3 .4  align=0. 8 
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material  material=GaN  arichn=24  arichp=96  edb=0.025  eab=0. 160 
material  material=GaN  permittivity=9.5  vsatn=2e7 
# 

model  material=AlGaN  conmob  fldmob  srh  print 
model  material=GaN  conmob  fldmob  srh  print 
# 

impact  bn=3.4e7  an=2.9e8 
impact  bn=3.4e7  an=2.9e8 
# 

contact  name=gate  workfun=5.0 

# 

#  SECTION  4:  Id-Vd  calculation 

# 

method  gummel  newton  itlim=20  trap  maxtrap=10  vsatmod.inc=0.01  carriers=2 
output  con.band  val.band  band.param  e. mobility  h. mobility 
solve  vgate=0 
# 

save  outf=NPSGaNHemtO_0.str 
tony  plot  NPSGaNHemt0_0.slr 
# 

log  outf=NPSGaNHemt0_0.1og  master 

# 

method  newton  trap  itlim=3  5  maxtrap=6  carriers=2 
solve  vdrain=0.0  vstep=0.1  name=drain  vfmal=  1  vsource=0.0  vgate=0.0 
solve  vdrain=  1.0  vstep=0.5  name=drain  vfmal=5  vsource=0.0  vgate=0.0 
solve  vdrain=5.0  vstep=0.75  name=drain  vfmal=20  vsource=0.0  vgate=0.0 
# 

save  outf=NPSGaNHemtO_Lstr 
tonyplot  NPSGaNHemt0_0.log 
# 

quit 
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